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Abstract. How sensitive are the predictions of inflation to pre-inflationary conditions when 
the number of efolds of inflation is not too large? In an attempt to address this qnestion, 
we consider a simple model where the inflationary era is preceded by an era dominated by a 
radiation flnid, which is conpled to the inflaton only gravitationally and which extends back 
to the Planck era. We show that there is a natnral generalized Bnnch-Davies vacnnm state 
for pertnrbations to the conpled inflaton-gravity-flnid system at early times. With this choice 
of initial state the model predicts interesting deviations from the standard power spectrnm 
of single held slow-roll inhation at large scales. However, the deviations are too small to be 
observable in near fntnre CMB observations. 
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1 Introduction and summary 

1.1 Background and motivation 

Inflationary cosmology as we know it today was formulated in the 1980’s [1~4]. This theory 
was designed to explain some of the key features of the Universe on large scales while reducing 
the amount of fine tuning required in initial conditions. In particular it provides a natural 


- 1 - 




explanation for the observed smallness of the Universe’s spatial cnrvatnre, the observed ho¬ 
mogeneity and isotropy on large scales, and the non-observation of magnetic monopoles [5-8]. 

The standard picture of inflationary cosmology, better known as single held slow-roll 
inhation, consists of a single scalar held (the inhaton) endowed with a canonical kinetic term 
and a slowly varying potential. The power spectrum of scalar perturbations predicted by this 
theory is largely consistent with the measured anisotropies of the cosmic microwave back¬ 
ground radiation (CMB) [9-11]. The measured temperature and polarization anisotropies are 
consistent with perturbations that are adiabatic, nearly Gaussian and nearly scale invariant, 
in excellent agreement with the predictions of the standard picture of inhation; see, e.g. Ref. 
[12] for a recent review. Additionally, the standard picture predicts a stochastic background of 
gravitational waves [13]. Primordial gravitational waves leave their imprint on CMB photons 
by inducing B-mode polarization patterns. Future detections of this effect could determine 
the energy scale of inhation. 

How robust the standard inhationary predictions are to effects of pre-inhationary phases, 
if any? If the currently observable modes exited the horizon many e-folds after the start of 
single held slow-roll inhation, then it is unlikely that any pre-inhationary phases can inhuence 
the dynamics of these modes. On the other hand, if the horizon exit of these modes occurred 
instead just after the start of inhation, is is conceivable that observational signatures of the 
pre-inhationary phase could be encoded in the longest lengthscale modes, those that exited 
the horizon hrst. 

There is a well-known discrepancy between the predicted and observed spectrum of CMB 
temperature anisotropies for multipoles Z < 40, which is small but statistically signihcant in 
the recent Planck data. Although speculative, it is conceivable that this discrepancy is due to 
evolutionary phases that preceded the single-held, slow roll phases. Several different scenarios 
of this type have already been explored in the literature, including preceding multiheld phases 
[14], initial non-slow-roll phases [15], and quantum mechanical tunneling from other vacua in 
a string theory landscape [16, 17]. Some of the effects observed in these scenarios were also 
found in single held models where perturbations are initially in non-vacuum states [18]. 

In this paper we address the question of the robustness of the standard inhationary 
predictions using a simple yet somewhat speculative toy model. We assume that the standard 
inhationary phase is preceded by an era of power law expansion driven by a perfect isentropic 
huid, whose equation of state is that of radiation The huid and inhaton are coupled only 
gravitationally. In Sec. 2 we discuss the properties of the background dynamics of this model, 
highlighting the existence of two distinguished transitions in the background evolution. Next, 
in Sec. 3 we derive the equations of motion for the scalar and tensor perturbations. In Sec. 
4 we discuss the adiabaticity of the perturbations in the distant past. We then argue for the 
suppression of the entropic perturbation modes in Sec. 5. We conclude in Sec. 6 and 7 with a 
few numerical case studies of the spectrum of tensor and scalar perturbations in this model. 

1.2 Summary of results 

• There are two important transitions in this model. The hrst occurs when the energy 
density of the huid coincides with the energy density of the inhaton, and signals a 
change from power law expansion to quasi-exponential expansion. The second, subse¬ 
quent, transition occurs when the breaking of the deSitter time-translation symmetry 
(parameterized by the time derivative H of the Hubble parameter) changes from being 

^This cosmological model has a long history. See for instance [19-21] for the study of the effects of this 
cosmological phase transition on the stability of quantum fields. 
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predominantly dne to the flnid to being predominantly dne to the inflaton. After the 
second transition the dynamics of the pertnrbation modes becomes indistingnishable 
from that of standard inflation. 

• Using the analytic backgronnd solntions fonnd in Sec. 2.3, we show in Sec. 4 that for 
snfhciently smooth inflationary potentials and ordinary flnids (p + Sp > 0), the flnid and 
inflaton pertnrbations deconple at early times. We also show that if the flnid has an 
eqnation of state parameter tc = 1/3 (radiation), all the pertnrbation modes of interest 
are in an adiabatic regime at early times, and that therefore one can dehne a generalized 
notion of the Bnnch-Davies vacnnm for the pertnrbations. 

• We show in Sec. 4 that the flnid and inflaton pertnrbation modes of interest deconple 
after the second transition. Therefore, both the backgronnd cosmology and the dy¬ 
namics of the pertnrbations become identical to single held slow-roll inhation after the 
second transition. 

• Onr model is incapable of generating signihcant entropic pertnrbations, an intnitive fact 
which we establish in Sec. 5. Thns effectively all pertnrbations are adiabatic. 

• We compnte the power spectra of scalar and tensor pertnrbations in the model, and 
hnd interesting deviations from standard predictions at long lengthscales. In particnlar 
there are oscillatory featnres in the low I end of both spectra (similar to those seen 
in trans-Planckian models), whose amplitnde depends on the ratio of the hnid energy 
density to the inhaton energy density at the epoch when the I = 2 mode left the horizon. 
In addition the scalar and tensor spectra are snppressed at low mnltipoles. 

• For some valnes of onr model parameters, the observed modes do not have a nniqne adi¬ 
abatic vacnnm at early times, becanse the mode evolntion has not yet become adiabatic 
when followed backwards in time before the Planck scale is reached. In this portion of 
parameter space, onr model loses predictability. By restricting the valnes of the model 
parameters, we can avoid this portion of parameter space. However, when we do so, 
the deviations from standard predictions are constrained to be small, of order a few 
percent. 

• We have not performed an analysis of the detectability of onr predicted effects. However, 
rongh estimates snggest that in the regime where onr model is predictive, the effects 
are too small to be seen in near-fntnre CMB experiments. 


2 The fluid-inflaton model 


We consider a cosmological model for before and dnring inflation consisting of a perfect 
isentropic flnid and a scalar inflaton held ip that are conpled only gravitationally. The action 
for the model in the metric signatnre (—h -|—1-) is 


S 






( 2 . 1 ) 


where R is the Ricci scalar, £/ is the hnid Lagrangian, V{'~p) is the inhaton potential, and 
we are nsing nnits with h = c = 1 = 8 ttG. The corresponding eqnations of motion are the 
Klein-Gordon eqnation 

VaV“(^ + V^{(p) = 0, 


- 3 - 


( 2 . 2 ) 




the Einstein equation 


(2.3) 


/o _ ^scalar i n-\ 

^ab — J-ab “T 


' ab 


where 

- ^gab[{V^? + 2V{^)\ 

and 

= {pf +Pf)UaUi,+Pfgab 

are the inflaton and fluid stress-energy tensors, and the fluid equation 

y^j^fluidab ^ 

Here Pf,Pf and are the fluid density, pressure and 4-velocity. We will later specialize the 
equation of state pf = Pf{pf) of the fluid to be that of radiation, pf = /0//3. 

In this section we will analyze the background cosmological solutions of this model, and 
in Sec. 3 below we will analyze the behavior of perturbations. 


(2.4) 

(2.5) 

( 2 . 6 ) 


2.1 Background cosmology 

For a homogeneous, isotropic and spatially flat background spacetime we write the metric as 

ds^ = —a{vi)^[dg^ — (dx)^], (2.7) 

where a{g) is the scale factor and g is the conformal time. We write the background solutions 
for the inflaton and the fluid as p{g) and p^{g), with the background fluid pressure being 
P/ = PfiP^f)- These solutions obey the two Friedmann equations 

+ r(i^) (2.8) 


4 {ml + w") = -p) - + F(^). (2.9) 

where primes denotes differentiation with respect to g and PL = a'/a is the conformal Hubble 
parameter, and the Klein-Gordon and fluid equations 


(f" ‘IT-Lp' -|- a^H^((^) — 0, 


( 2 . 10 ) 


p'f + 3PL{p^f TT*/) ~ O' 

The continuity equation (2.11) yields 


0 ^ ^0 
W 03(1+22;) ’ 


( 2 . 11 ) 


( 2 . 12 ) 


for an equation of state pf = wpf, and a constant of integration Sq. In our analysis of the 
background dynamics we will mainly use the first Friedmann equation (2.8) and the two 
equations (2.10) and (2.11). 
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2.2 The background dynamics and parameters 

The background fluid-inflaton Lagrangian depends on a number of parameters, namely the 
equation of state parameter w of the fluid and all the coupling parameters in the inflaton 
potential. In addition, the space of solutions of the differential equations (2.8), (2.10) and 
(2.11) is four dimensional, requiring four initial conditions to determine a solution. One of 
these parameters is gauge, since we can resale the coordinates to set the value of a at the 
initial time to unity. The remaining three parameters can be taken to be the initial value 
and initial time derivative of the inflaton, and the initial ratio PipjPf of the energy densities 
of the inflaton and fluid. 

In exploring this space of model parameters, our strategy will be to choose models that 
agree with observations in the absence of a fluid, and then gradually modify the model by 
dialing the parameters related to the inflaton at some initial time when the energy density 
of the fluid is Planckian. In particular we will restrict attention to models in which, in the 
absence of a fluid, inflation occurs at the GUT scale of ~ 10^® GeV. We will use as an example 
the large field mode V = vn^ip^ j2 to construct numerical examples of the power spectra of 
the scalar and tensor perturbations in sections 6.1 and 7.1. 

To illustrate the qualitative features of the the background dynamics. Fig. 1 shows the 
evolution of the fluid and inflaton energy densities, the Hubble parameter H = T-L/a, and the 
inflaton ip. For this example we have taken V{p) = w?p^I2 with m = 10“^, and for the fluid 
we have taken tc = 1/3 (radiation) and used a Planckian initial energy density. The initial 
inflaton kinetic and potential energies are 

^ .. V{p) ^ 10-«, (2.13) 

which is about the grand unihcation energy density (note that as stated previously we are 
using units with nip = 1). As seen in Fig. 1, the background dynamics starts in a radiation 
dominated era. The energy densities become equal at N ~ 5, after which the energy density 
of the radiation keeps decreasing rapidly while the background dynamics becomes similar to 
that of single held slow-roll inhation. The ensuing near constancy of the Hubble parameter 
bears evidence to the Universe expanding quasi-exponentially, as expected from the standard 
picture of inhation. 



Figure 1: Left: Logarithm of energy density p as a function of number N of efolds of inflation. The radiation 
fluid is shown in black and the inflaton in red. Note that the radiation fluid is initially dominant but rapidly 
becomes subdominant. Middle: Log of the Hubble parameter as a function N. Once the inflaton starts 
to dominate at ~ 5, the Hubble parameter becomes approximately constant and the expansion becomes 
quasi-exponential. Right: The inflaton field p as a function of N. Note that the evolution of the inflaton is 
nearly frozen during the fluid dominated era. 

During the early, huid dominated era, the inhaton held evolves very little and is effec¬ 
tively frozen. This can be seen in the third panel of Fig. 1. The freezing of the inhaton is 
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easy to understand: defining N = Ina, we have dlmp/dN = (p/{(pH). When pj ^ p^, the 
fluid dominates the right hand side of the Friedman equation (2.8), giving 

dlnp p ^ 

~dir^ 

Thus we have dlnp/dN <C 1 whenever p^ <C Pf, as long as the inflaton held is larger than 
the Planck scale. One can also think of this effect in terms of an enhanced Hubble damping: 
the Hubble damping term in the Klein Gordon equation (2.10) is much larger than it would 
be without the huid, since is much larger than p^. Because of this enhanced damping, the 
inhaton kinetic energy rapidly becomes much smaller than its potential energy. For example, 
in the model shown in Fig. 1, the kinetic and potential energies are initially equal, but by the 
time of huid-inhaton equality the kinetic energy has decreased by more than three orders of 
magnitude. This then implies that the inflaton dominated phase starts out in the slow-roll 
regime: the value of the slow roll parameter = p ''^/when the inflaton begins to 
dominate is 10-3. 



Figure 2: The inflaton slow roll parameter as a function of number N of efolds, for the fluid-inflaton 
model with potential V{p) = rrdp^j2 (blue), and for the corresponding inflaton-only mode (dashed red). 
In this example, the initial fluid energy density is the Planck scale, while that of the inflaton is about the 
grand unification scale 10“^. The slow roll parameter in the inflation-only model starts at some value set by 
the initial conditions, and rapidly decays to the asymptotic track determined by the potential. By contrast, 
in the fluid-inflation model, the slow roll parameter is initially driven to a very small value during the fluid 
domination regime, and then relaxes upwards to the asymptotic track determined by the potential in the 
inflaton dominated regime. 

As discussed in the introduction, there are two important transitions that occur in the 
ffuid-inffaton model. The hrst is the transition already discussed of p^^ ~ p/, where the 
background expansion changes from the fluid dominated power law expansion to the inflaton 
dominated quasi-exponential expansion. After this hrst transition there is an approximate 
de-Sitter expansion. 

At the hrst transition, the comoving horizon scale Lc = 'HT^ has a local maximum. It 
increases during the huid dominated era, since 

-^Ccoc p + 3p, (2.15) 


(2.14) 

P\l Pf 
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Figure 3: The comoving scale aH versus number N of efolds for a fluid with an equation of state parameter 
w = 1/3 and an infiaton with a potential. The two important transitions are shown. 

aH 



and during this era p and p are dominated by the fluid which we assume obeys the strong 
energy condition w > —1/3. It decreases during the inflaton dominated era as usual. Because 
of the existence of a local maximum, there exists a minimum comoving wavenumber /Cmin such 
that modes with k > k^am leave the Hubble horizon during the inflaton dominated phase, but 
modes with k < kmin are outside of the Hubble horizon throughout the fluid and inflaton 
dominated phases. Using 

= (2.16) 

(where we have approximated ~ Aq at early times), minimizing with respect to a and 
noting that we choose i5o = 1 at a = 1 gives 

l+3m 

fcmin ~ (2.17) 

For example, for radiation with w = 1/3 we have fcmin ~ 

The comoving Hubble horizon defined above is completely determined by the background 
dynamics. However, in this model there are different scales analogous to the Hubble horizon 
that are more relevant in discussing the evolution of the perturbations. To see this, we use 
some of the results derived in Sec. 3 below. Let us initially neglect the interaction terms 
between the fluid and the inflaton perturbation modes in Eq. (3.27). We see that each mode 
obeys a harmonic oscillator equation with the effective masses m^(fe, rj) and m‘j{k, rj) given by 
the diagonal elements of the matrix f2 defined in Eq. (4.2). One can see that the perturbations 
evolve differently in two distinct regimes. The first regime is when ~ ~ k'^, and the 

second regime is when k'^ is much smaller than |m^ — k‘^\ and |m^ — c^k‘^\- The transition 
occurs roughly when 

k'^ ~ |m^ — or — k'^\. (2.18) 
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Based on this transition property, we define 


£ 


c—inflaton 


1 

V\k^-nu 


and 


C 


c—fluid 


- n22/cj 


(2.19) 


in analogy with the Hubble horizon. 

Nevertheless, we will indicate in Sec. 4.1 that the scalar modes with wavenumbers larger 
than /cmin are the ones that are of primary interest to us. These scalar perturbation modes are 
initially in an adiabatic regime. Furthermore, it can be shown that after the first transition, 
while the background geometry is nearly de Sitter, the above defined adiabatic horizon scales 
are within a factor unity of the Hubble horizon. Thus anytime we refer to a mode exiting the 
horizon we shall simply mean its wavenumber k becoming equal to T-L. See Fig. 4 for a plot 
of all three horizon scales. 


Horizon ^ 



Figure 4: The inverses of the three horizon scales Lc, dc-inflaton and £c-fluid versus number of efolds N. The 
green plot is the plot of = T-L, which is a background function used to define a local-in-time horizon. The 
dotted blue plot is the plot of the inverse of the inflaton adiabatic horizon scale , and the dashed 

black plot is the plot of the inverse of the fluid adiabatic horizon scale d“2fluid associated with a fluid with an 
equation of state parameter ui = 0. The adiabatic horizon scales roughly signify the boundary between the two 
distinct regimes for the perturbations; namely the highly oscillatory regime and usually the non-oscillatory 
regime. We have also included a plot of the inverse of the fluid adiabatic horizon scale for the case of radiation 
w = 1/3 (the red plot vaguely seen at the bottom). Note that the radiation case is an exception due to being 
conformally invariant in any (conformal) power law cosmology. 


We now turn to a discussion of the second important transition in the model. This 
transition occurs when the breaking of the time translation symmetry (characterizing the 
deviation of the spacetime from the de Sitter spacetime) switches from being predominantly 
due to the fluid to predominantly due to the inflaton. Recall that 



a 



(1 + w) 


(j3(l+u;) 


o2 J 


( 2 . 20 ) 


where H is the Hubble parameter and dot denotes differentiation with respect to the coor¬ 
dinate time variable t. Before the transition, the first term in Eq. (2.20) dominates, and 
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afterwards the second term dominates . One expects the background dynamics to begin to 
converge to that of single held inhation once the time translation symmetry is broken by the 
inhaton kinetic energy. We can also quantify this transition in terms of the huid and the 
inhaton contributions Sf and to the primary slow roll parameter e = given by 

( 2 . 21 ) 

The second transition occurs when £p ^ Sf. To be more precise, we note that slow-roll 
inhation requires the smallness of another parameter known as the second slow-roll parameter 
r] = —H / {2HH). This parameter ensures the smallness of the inhaton held acceleration. Like 
the primary slow-roll parameter e, this parameter also has contributions due to the huid, r]f, 
as well as the inhaton, ijp, that are given by 




Pf 

eiL2’ 




eH^ 


( 2 . 22 ) 


We expect to recover a single held slow-roll background dynamics by the time e ~ and 
r] ~ 7]p. Thus, we take these two conditions to be the dehning conditions for the second 
transition. These conditions together place an upper bound on the energy density of the huid 
at the moment of the second transition, that is 


Pf ^ EpTipH . (2.23) 

In what follows we shall take Pf = SpTjpH'^ as the dehning moment for the second transition. 
The huid energy density falls below the upper bound of Eq. (2.23) at an e-fold around 


log 


3(l+m) 




(2.24) 


where Ai is the approximate value of the inhaton energy density during inhation. For a huid 
dominated phase with w = 1/3 and an Ep = \rip\ = 0.01, we hnd that after N ~ 6.9 e-folds 
the background dynamics is practically identical to that of a single held slow-roll inhation. 


2.3 Analytic solutions to the background equations in asymptotic limits 

The background Eqs. (2.8) to (2.10) can be solved analytically in the limits when either 
one of the energy densities is sufficiently subdominant to the other one. The strong Hubble 
damping imposed on the motion of the inhaton held up until the second transition allows one 
to approximate the energy density of the inhaton as a cosmological constant to leading order. 
This means that in either regime, we expect to approximate the inhaton potential using the 
hrst few terms in its Taylor series expansion 

V{(p) = A* -h V'{ip^){ip - ¥?*) + 0[{ip- ¥?*)^], (2.25) 

where is the value of the inhaton held either at the beginning of the huid dominated era, 
or at the onset of the second transition, when the inhaton energy density becomes dominant. 
This property will be particularly useful when we discuss the decoupling and adiabaticity of 
the perturbation modes at early times in Sec. 4. 

Based on this insight, we regard the background cosmology at early times to be the huid 
plus a cosmological constant and we study the perturbations around this model induced by 


9 










the variations in the inflaton field. Examining Eqs. (2.8) and (2.10) in this limit reveals that 
we can solve the background equations perturbatively using the following expansions for the 
scale factor and the inflaton field 


a = a [l + + 0{a ‘^)\, (^ = (/pq [l + + 0((T^)] , (2.26) 


where tpQ is the initial value of the inflaton field, a is the scale factor due to the dominant fluid, 
and a is the perturbation variable taken to be the dimensionless ratio of the energy density 
of the inflaton over the fluid when the fluid is Planckian. It turns out that the background 
equations are easily solved using this perturbation scheme. What makes these equations 
particularly simple to solve is the fact that to each order in cr, the Eqs. (2.8) and (2.10) 
decouple. In other words, the order term in a is used to solve for the fi” order term in 

Lp. The reason is simple to see; as we commented earlier, this perturbation method is really a 
perturbation around the fluid plus cosmological constant model generated by the variations 
in the inflaton field. Indeed the first order correction to the scale factor, is due to the 

zeroth order inflaton field acting as a cosmological constant with p^p = E(y?o) = ^o- For fluids 
with w > —1/3, one finds that decays like ryi+a™ as r] (or a) approaches zero. 

Using the zeroth order scale factor a, one then finds that there exists a solution for which 

/I \ 6+6w 

decays like ryi+s™ . This confirms that the expansion (2.25) is valid at early times. Other 
solutions for diverge in the limit of p going to zero, spoiling the fluid domination at 

early times. We therefore specialize to this class of solutions for the inflaton field. 

To provide a concrete example, we apply the above formalism to the case of a radiation 
fluid [w = 1/3) coupled to an inflaton with V{(p) = rr?tpP‘/‘I. In this case, we find 


a{j]) 


1 + 


4 iplihpl - 12)m^Tg 8 , 
180 ' 194400 ' 


v>{v) = v>0 


1 - 


(2^o-15)m%^ 8 I 
194400 ' 


(2.27) 


Next, we analyze the background cosmology after the second transition, when the back¬ 
ground dynamics begins to converge to single field slow-roll inflation. As a result of the 
initially-fluid-induced Hubble damping, the inflaton field behaves as a cosmological constant 
to a very good approximation for a while even after the second transition. Thus, we can em¬ 
ploy a similar perturbation scheme to study the perturbations to this cosmological constant 
model induced by both the inflaton field variations and the remaining fluid energy density. 
To do this, consider the following expansions for the scale factor and the inflaton field, 

a = a[l-|-e(^AaF) +0{e^)], p = ipi +0{e^)], (2.28) 

where pi is the value of the inflaton field at the second transition, a is the scale factor due 
to the dominant cosmological constant, and the perturbation parameter e^, is the inflaton 
slow-roll parameter introduced in Eq. (2.21) that is evaluated at the second transition. As 
mentioned in the last subsection, this parameter breaks the time-translation symmetry during 
the ensuing slow-roll inflation. 

To proceed, we have to determine at what orders in the added fluid will contribute 
to the perturbative expansions in Eq. (2.28). This can be understood using the upper bound 
derived for pj in Eq. (2.23). There we found that once pf ~ e^Ai, the background dynamics 
begins to converge to that of single field inflation (we have taken for simplicity). 
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This regime is precisely where our perturbation scheme is valid. It can then be seen from 
Eqs. (2.8) and (2.10) that the leading fluid corrections to the scale factor and the scalar field 
are at most of order e^. 

It actually turns out that the leading fluid correction to the scale factor is of order 
while its leading correction to the scalar field is of order e^. This is a result of the same 
decoupling of Eqs. (2.8) and (2.10) to each order in perturbation that also occurred in the 
distant past regime. In this case, the decoupling occurs because the zeroth-order inflaton field 
(the cosmological constant) induces variations in itself via the Klein-Gordon equation (2.10), 
i.e. the first order scalar field correction is generated by the cosmological constant 

V{ipi) = Ai. Therefore, it is easy to see that the order term in a can be used to solve 
for the ^ order term in ip. 

For the case of a radiation fluid (w = 1/3) coupled to an inflaton with V{ip) = /2 

a straightforward calculation gives 



«(f?) = \l ^ 


1 - 


27 — 2 IA 7 / — TATy"^ + Ary’^ + 42Ary log | A?y| 


42A7y 


£in + 


(p(r]) = pi 


1 + 


2 — 2Ary^ + 6 log |A7y| _ 


6 


£10 + 


(2.29) 


Here Ai = rr?p\l2 and Ary = ry — r/end) with ryend being the conformal time at which the scale 
factor a = {\fKi/\r]) becomes infinite. Also, in the above calculation we have set Ary 

equal to one at the second transition. 

We will use the equations derived here in Sec. 4 to show that the radiation and inflaton 
scalar perturbations will dynamically decouple in the distant past (radiation dominated) and 
the future (inflaton dominated) regimes. 


2.4 Post-inflation dynamics: reheating 

Conventionally, inflation ends when the primary slow-roll parameter e becomes equal to unity. 
At this time, the Universe no longer expands quasi-exponentially since the kinetic energy and 
the potential fluctuations of the inflaton field are large enough to completely break the de 
Sitter symmetry. Once this occurs, the Universe enters a phase commonly known as reheating. 
During this phase, the inflaton field decays to other particles, though the exact mechanism 
remains unknown and is largely unconstrained [See [22] for a good review]. 

The simplest approach to understanding reheating is to assume that inflaton decays to 
a single component fluid with an equation of state parameter Wre- The decay occurs while 
inflaton, now in the form of non-relativistic particles, oscillates at the bottom of its potential. 
The Klein-Gordon equation (2.10) is then corrected to take into account the decay process, 

p'^ + ^T-Lpip = -aVpp, (2.30) 

where T is the decay width associated with an inflaton decay process. 

In this work our interest in reheating lies in the contribution of the reheating phase to 
the overall redshift at which the energy density of the fluid pf becomes Planckian. This allows 
us to determine the range of physical scales at the initial time that corresponds to the range 
of GMB observable scales today. It turns out that we can write the ratio of the initial scale 
factor Oj to today’s scale factor oq as 

Oi g—Mot ®end 

Q/Q Qj'p ^^0 


(2.31) 
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The above A^tot is the total number of e-folds between pf = 1 to the end of inflation, Oend 
is the scale factor at the end of inflation, and is the scale factor at the beginning of the 
radiation dominated era. It is a straightforward computation to show that 

®end _ ( \ 3(1+Wre) 

^ Pend ^ 

where pr is the energy density of the radiation and pend is the total energy density at the end 
of inflation. Similarly, the ratio a^/ao is computed with the cosmological parameters provided 
by the flat ACDM model. 

Here we shall restrict attention to the most trivial theory of reheating, i.e. the in¬ 
stantaneous reheating, when constructing numerical examples in sections 6.1 and 7.1. For 
instantaneous reheating Oend/or = 1- However a better motivated model of reheating, such 
as the canonical reheating, has Oend/or different than unity. In the canonical reheating model 
the inflaton decays to a single component non-relativistic matter. The duration of reheating 
in this case is determined by the energy scale at the end of reheating. This energy scale 
is proportional to the so-called reheating temperature Tre, which is the temperature of the 
radiation fluid at the beginning of the radiation dominated era. When taking into account the 
constraints placed on inflationary models coming from reheating considerations, one needs to 
consider such effects for various Tre [See [23, 24] for a recent study on this subject]. Although 
as it turns out, a Tre lower than the energy scale at the end of inflation (as in the canonical 
reheating case) only renders our predictions less interesting. 


(2.32) 


3 Dynamics of the perturbations 
3.1 Fluid action 

In this section we derive an action principle for the coupled fluid-inflaton-gravity perturba¬ 
tions. We start by describing in more detail our action principle for the fluid. 

We describe the fluid flow by a set of comoving coordinates x“(a*,A). Here a* are a 
set of Lagrange coordinates chosen on an arbitrary spacelike hypersurface and A is an affine 
parameter, which together provide sufficient data to label a fluid flow line passing through a 
given point in spacetime. The flow of a perfect fluid is constrained by the conservation of 
the number density current along all flow lines. 


Va = 0. 


(3.1) 


The above J“ is the number density current of the fluid dehned as nM“, with being the 
fluid four-velocity and n the fluid number density. One then defines the ffuid four-velocity by 
dx°‘/dX {—dx^/dX dxa/dX)~^/‘^, thereby arriving at [8] 


T{a^) dx 


where is some function of the Lagrange coordinates and J = T>{x°‘)/T>{a^, X) is the 

Jacobian associated with the coordinate transformation from the Lagrangian coordinates to 
the comoving coordinates. It can be easily seen that J dehned in Eq. (3.2) solves the 
continuity equation (3.1). 
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With this description of the fluid flow, the general action for this model in the metric 
signature (—h ++) becomes 


S 


d^x 


-R 

.'2 


+ ^r. 


(3.3) 


where g is the determinant of the spacetime metric, R is the Ricci scalar, and Cm is the 
matter Lagrangian density. For convenience, we have set h = c = mp = 1, where nip is the 

reduced Planck mass dehned as (SvrG)”^/^. The matter Lagrangian density for this model is 
2 

Cm[gab,C,x'^] = -pf{\J\) - -y{c)- (3.4) 

The above p/ is the energy density of the fluid and (p is the inflaton held with a potential 
V{ip). Also (V(^)^ is a shorthand for gah^°'p^^<p as usual. Note that the huid pressure pf 
can be dehned in terms of the huid energy density pf via [25] 

Alternatively, a huid with an equation of state parameter w has its energy density and pressure 
related hy pf = wpf. 


3.2 Scalar perturbations 

In this section we will derive the linear equations of motion for the gauge invariant variables 
associated with the scalar perturbations to the background huid-inhaton model introduced 
in Sec. 2. It is convenient to derive these equations by expanding the action (3.3) to second 
order in the scalar perturbations. Working with the action has privileges over the perturbed 
linearized Einstein equations. The clarity of the process of integrating out the non-dynamical 
components of the metric as well as the quantization of the perturbations are the main 
advantages of working with the action principle. 

The scalar perturbations to the FRW metric in an arbitrary gauge can be written as [7] 


ds^ 


-aiv? 


(1 + 2(l))drf‘ 


2B^idx'^drj — [(1 — 2ip)5ij + 2E^ij\dx^dx^ 


(3.6) 


for some scalar function </>, ip, B, and E. Similarly we write the perturbed scalar held as 


V?(r/,x*) = (fig) +Sp{g,x^). 


(3.7) 


We now insert the expansions (3.6) and (3.7) of the metric and inhaton into the action (3.3), 
and expand to second order in the perturbations. The hrst order terms must vanish identically 

^One formulation of the Lagrangian density of a perfect isentropic fluid is given by 

-V^Pfi\j\)-xVar, 

where A is a Lagrange multiplier. If we use the Lagrangian fluid variables explained on page 12, then = 0 

by virtue of Eq. (3.1). In this case, the fluid Lagrangian density reduces to —^—gpf (|7|). We refer the reader 
to Sec. 5 of [25] for an extensive discussion of the Lagrangian fluid variables. The reader might also find various 
common formulations of the fluid action principle discussed on pages 31-36 of [25] insightful. 
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when the background variables are on-shell, so we drop all the hrst order terms. The second 
order gravitational and inflaton terms are 


62 iSgravity 


= j dr] j + + 

-An{(p -k - E') + m-ilr'V‘^E - - E') - mV^lr-VB 

+ 'ip)V^E - 4nV'^EV‘^{B - E') + AnV'^EV^B + m'^iV'^Ef 

+?>n^{VBfY 


V^|;) 


(3.8) 


and 


*525*inflaton 

= j dv j d^x^{{a‘^Viip)-^){^^-3^lj‘^-{VBf + (y‘^Ef + 2^ljV^E + 6(J)i; 

-2(j)V^E) - 2(((. - SV’ + V^E) {a\^{ip) 6ip + (^'[V - 6ip']) - [(Vht/?)^ 

+2ip'VB.V6^ + + Ac^ip'dip' - }, 

(3.9) 

where = S^Widj and .VB = d^^A^iBj. 

Perturbing the fluid Lagrangian is a more delicate task. Note that the energy density of 
the fluid pf itself is not a dynamical variable. For a perfect isentropic fluid, the energy density 
is a function of the fluid number density re. One needs to perturb the fluid Lagrangian by 
perturbing re, and use (3.1) to relate the perturbations in re to the metric perturbations and 
the perturbations in the fluid four-velocity, the latter being related to the perturbations of 
the comoving frame of the fluid. Following [8], we denote the perturbations to the comoving 
frame of the fluid by a shift vector ^*( 77 , x-^). Since we are interested in scalar perturbations 
here, we define = di^ for some function In short, we are writing the number density 
perturbations, 5n{ri,x^), in terms of ^(r/, x*) and other metric perturbation functions. The 
details of this calculation can be found in [8]. The final result is 


<525'fluid 




— 3i^?/i + 4>V‘^E 


^PV^E+'^{S7^Ef - E^^jE^ij 


+ \{VBf'^ + {pj + pj) + B4[, + 

-Ic2(p0 +p0)(37/> - V^E - V2C)'|, (3.10) 


where Cg is the sound speed of the fluid defined by Adding the contributions 

(3.8) to (3.10) and using the background equations of motion (2.8) and (2.9) together with 
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integrations by parts in time and space, the total variation in the action simplifies to 


628 = JdT] j y{+ + + 

+ 2 (VV ')2 + a^{p) + B+ a^(p‘} + pj)(2(/>V^E - 6</)^ + 2c/>V^^) 

-cla^{p) +p|^)(3V' - + ^.'202 + _ 3 ^ + V‘^E){p'5pa^)' - 

-2p'VB.V6p - - ^p'5p' + (V)^}- (3-11) 

In Eq. (3.11) six pertnrbation functions appear; fonr associated with the metric, one associ¬ 
ated with the scalar field, and one associated with the flnid. It is convenient to replace the 
flnid pertnrbation fnnction ^ with a three-velocity potential defined as [8] 

e=‘^^{e + B) (3.12) 

Cs 

which is related to the pertnrbation of the fonr-velocity 6u^ of the flnid by 

9i = -‘^^^5ui. (3.13) 

’ Cs 

In Eq. (3.12) /3 is defined as (a'^/2)(p^ +Pj)- Using the definition (3.12) and the constraint 

- V^E - V2| + 4) (3.14) 

Z Cg 

obtained from varying the action (3.11) with respect to we eliminate all and ap¬ 

pearing in Eq. (3.11) in favor of 9 and the other metric pertnrbation fnnctions. The final 
resnlt is 


628 = J dp J d^x + n 4 )f+ 2 ^(j)^- 2 a^Vi^.{ 2 V 4 >-Vilj) 

- lcl{V9f - Qcsa^9i^' - 2a^Cs9l^^ - Aa^V\B - E') 

+ Bcf) - E + a2(^'2())2 + 2{4> - ^i>){(p'6(pa^y - a^{V6(p)‘^ 

—a‘^V^y,S(p'^ — Aa^cjiip'5(p' + a'^{Sp')‘^ . (3.15) 

Note that the action (3.15) depends on the pertnrbation fnnctions B and E only throngh the 
combination B — E'. Varying the action (3.15) with respect to B — E' gives 

'll;'+ n(j) =\ip'5p - ^y/^9. (3.16) 

2 2a 

As can be seen from Eq. (3.15), 4> is non-dynamical and it can be integrated ont nsing Eq. 
(3.16), i.e. (j) can be eliminated in favor of ip, 9, and dip. Observe that by doing this we 
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eliminate E, B, and (j). This way we are left with three dynamical fields; V’j and ^4^- The 
action (3.15) after integrating out B, E, and cj) becomes 


52 S =[ dr] [ d^x- 


3 2o {CsV^o)'"^ {CsV^o)' 


— —cB + —— — - -h 

2^^ 2n^ 2%'^ 2'H‘^ 

, ( 2a/3^/^ 


+ 


c^jda? 


2csy/^a 
/2 


, Hsf W 2a\ , fl,//2a/3T2 2(cs^/^a)' aif 


Uci 

'2a^Vu 


2\/^Ofl/ ,/ , ,/2/2a2/3 oV'^\ , lfl/2 , 2,i: /^2 , , 2^v7/^2 

(^ + ^) + 2^ + — VV^.VV^ + 2a (VV^) 

+2:^c, - ^cliVef - a\V5^f + V 




+ 


2 /2 \ 
~m 


(3.17) 


The functions 5(p^ 4;, and 0 form the following two gauge-invariant variables 


u = ^(6» - 2z4}), 


vr = a5(p + C'4’, 


(3.18) 


where z = a^/~0/{'Hcs) and ^ = aip'/Ti. The function u is the gauge invariant fluid velocity 
potential, and the function vr is the gauge invariant inflaton perturbation. 

So far in this computation we have not made any gauge specialization but kept the 
gauge arbitrary. One might imagine that this computation would be easier if one imposes 
some gauge conditions at the start, for example the Newtonian gauge conditions B = E = 0. 
However, this gauge specialization is too restrictive and does not allow us to properly integrate 
out the non-dynamical variables. With that gauge choice, we would not have been able to 
derive the constraint equation (3.16) needed to integrate out cj). Similar difficulties arise in 
the synchronous gauge where cj) = B = 0 (but not in the inflaton comoving gauge where 
Sip = E = 0). 

Rewriting the action (3.17) in terms of 12 and vr defined in Eq. (3.18) gives ^ ^ 


S 2 S / dr] j d X "^Tfiuicj -|- Tjnflaton T Tint ^; 


where 


Tfiuid = ^ ( - c2(Vz.)2 + 


z" 

-h 4H 

2 ; 


is the fluid Lagrangian density. 


(3.19) 


(3.20) 


Tinflaton 


tt'^ — (Vtt)^ -|- 


u c 


+ ^(4z^A-(zBy) 


vr 


(3.21) 


^The action presented in this form is not applicable for the case of pressureless dust fluid. A redefinition 
of variables is necessary for this special case. 

^For linearized scalar cosmological perturbations, isentropic fluids admit a simple scalar field description. 
For the case of radiation fluid, one can replace the fluid Lagrangian density —p/ with (V't/i)^/4. Cosmologies 
of scalar fields with non-canonical kinetic terms of this form had been previously studied (see e.g. [26]). 
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is the inflaton Lagrangian density, and 


C 


int 


= (4^/2A j - ^/2B^ 


rB , 
i^vr + V 2 —i^vr 


is the interaction Lagrangian density. Here 


A = 


,2 


!2 


sn 


8W2' •’ 2c, 2H 


(3.22) 


(3.23) 


and 


B = 


V' cj 2n^ 


(3.24) 


are functions of the background variables. We remark here that the fluid influences the 
inflaton perturbations in three steps: hrst through the explicit interaction terms given in Eq. 
(3.22), second by modifying the expressions for the background functions A and B that involve 
the fluid variables, and third by modifying the background solutions which also influence the 
inflaton perturbations. Finally, observe that turning off the inflaton by setting vr = 0 gives the 
standard gauge invariant Lagrangian density for the isentropic fluid [8], while turning off the 
fluid by setting ly = 0 gives the standard gauge invariant action for the inflaton perturbations 
[27-29]. 

We can write the action (3.19) more compactly as 


62S = J dri J d^x - 

where the matrices are dehned as 

n = 


n'^n' - c^v*n + rn + n'^ An 


rT„2v7i 


/T' 




1 0 
0 c) 


2 > 




0 


-C2f 0 



^ ^ ^+4A 


^4^/2A| - ^Bf - 


(3.25) 


(3.26) 


We have suppressed the dependences on rj for all the background functions above. The Euler- 
Lagrange equation of motion for this action is 

n" - [c^v^ + r - ^A'] n + An' = o. (3.27) 

Dehning the Fourier modes of the perturbation functions using 

j (3.28) 

the corresponding equations of motion for each mode become 


u" + - 4a')l' — (4\/2A- — V2B^')Tr — VI—'k' = 0, 


H I 

vr + 


V C 

^ - ^{4z‘^A-{zB)') 




-B 


c 


(3.29a) 


B' . n, .Z\ n,B 


c 


c 


vr + ( \/2— - 4 V 2 A- ) u + \/2—u' = 0. (3.29b) 


c 


We have omitted u and tt dependences on r] and k. This is a system of coupled oscillators 
with time dependent coefficients. 
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3.3 Tensor perturbations 

To derive the linear equations of motion for the gravitational waves produced in this model, we 
expand the first term of the action (3.3) to second order in the tensorial metric perturbations. 
Unlike the scalar perturbations, the metric tensorial perturbations do not couple to the matter 
perturbations. The tensor perturbations to the metric are defined by [5] 


ds‘^ = —a{r])'^ [dr/^ “ (1 + hij)dx'^dx^] , 


(3.30) 


where hij is a traceless and transverse matrix, i.e. = dihij = 0. It is not difficult to see 

that the tensor perturbations defined as such are gauge invariant [5]. We can express hij as 


hij{x,r]) 


S / 


d?k 

(27r)^/2 


a{v) 


efj(k) + c.c. 


(3.31) 


where A denotes the two independent polarization tensors which are traceless and trans¬ 
verse and satisfy = 2(5"^'®. 

The well known tensor perturbation action is given by [7] 


d25te 


1 

8 


drj 


dx^o? 


Kj - i^hijY 


(3.32) 


The equation of motion for the gravitational waves in the Fourier momentum space is given 

by 

h''^{k, rj) -|- u}‘^{k, r])h^{k, rf) = 0, (3.33) 

with uj‘^{k) = k'^ — a"/o,. This equation is shared in all models where the matter content 
is minimally coupled to gravity and the gravitational field action is the Einstein-Hilbert 
term. What renders the evolution of the gravitational waves different in various models is 
the difference in the background dynamics, i.e. the difference in a"/a ratio. From the two 
Friedman equations (2.8) and (2.9) we have 


^ = y (^ - 3p) = y [{Pv’ - ^P<p) + P/(l - 


(3.34) 


For the case of radiation {w = 1/3) coupled to inflaton, the second term in the bracket 
vanishes and a”/a is given by the same expression as in the case of the inflaton-only model. 
However, the background geometry and the dynamics of the inflaton field is still influenced 
by the presence of the radiation fluid. We will explore these differences in Sec. 7.1 where we 
numerically compute the spectrum of the primordial gravitational waves in the fluid-inflaton 
model. 


4 Past and future decoupling of the scalar perturbations 

An interesting property of the fluid and inflaton scalar perturbations in this model is that 
they are only effectively coupled while the background energy densities are within a few orders 
of magnitude of one another. To establish this, we must examine the equations of motion 
(3.27) both at early times and at later times when the background dynamics converges to 
that of single field inflation. The equations of motion (3.27) are of the form (after Fourier 
transforming to the momentum space) 

Hj -I- [Hjj -I- -Xij]Uj + XijWj = 0, (4.1) 
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Figure 5: Below is the plot of a" ja for the fluid-inflaton model (blue) and the inflaton-only model (dashed 
purple) with the if? potential. The initial values in the inflaton-only model are chosen such that the background 
dynamics becomes identical to the fluid-inflaton model after the second transition. 

a"/ a 



where 11 and A are defined in Eq. (3.26), and 


n = -clk^ - r = 


^ ^ [4z^A - {zB)'] -4\/2Af + 

\-4V2AI + ^B^ + ^f c^^k^-^-4A ) 


(4.2) 


As can be seen from Eqs. (4.1) and (4.2), Qn and Q 22 can be interpreted as the effective 
masses of the inflaton and fluid in this model. Turning off the interactions, they become 


rill —t — C”/( and Q 22 —t c^k'^ — z”jz 


(4.3) 


as they would be in inflaton-only and fluid-only models. From Eq. (4.1), one expects that 
the perturbations decouple once 

^12) ^12 + 2-^12riii,n22- (4.4) 

We can check this condition using the asymptotic background solutions obtained in Sec. 2.3. 
The result for the distant past limit is given in Table. 1. 



-T<W<k 

«;= 4 

rill, ri22 

cjk'^ + (constant) tj ^ 

cjk'^ + (constant) + -|- ... 

ni2 -|- 2-^12 (* j) 

7 8+6411 \ 

7 8 +(jm \ 

Olrj 1+3™ j 

0(771+3™ j 

^2 

0^771+3™ ^ 

0^771+3™^ 


Table 1: Distant past forms of the background functions appearing in Eq. (4.1). 


This calculation shows that the perturbation modes decouple in the limit r/ —?■ 0. Another 
question that could be asked is: does a mode of a given wavenumber k decouple before 
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the Planck scale is reached? 
dimensionless ratio 


To address this question, one could numerically compute the 


^12 + 5^12 


(4.5) 


just as pf becomes Planckian, and see whether it is small or large compared to unity. 

Next we would like to show that this dynamical decoupling of the scalar perturbations 
also occur when the background dynamics begins to look like single field slow-roll inflation, 
i.e. when pf < This can be done using the late time solutions of the background 

equations (2.29). The result is provided in Table. 2. From this we see that early on during 



ru = 4 

Hii - k'^ 


H22 - c^k'^ 

(4^-4Ar? + 2A7?4)e-^ + 0(4) 

U12 + \X'i2 


A'i2 



Table 2: Future forms of the background functions appearing in Eq. (4.1). 


the second transition we have 

rill — ~ ~ 1 > ^12 + -A'i2 ~ A 12 ~ 

r!22 - ~ 0(4) < ^12 + ^a;2 ~ a;2 ~ (4.6) 

where we have used Ap = 1 at the second transition consistent with our choice of scale factor 
normalization in Sec. 2.3. From Eq. (4.6) we see that the inflaton mode has effectively 
decoupled from the fluid mode by the second transition. However, this does not appear to be 
the case for the fluid mode. It is in fact simple to show analytically that 

z" I 3 1 

^22 — =-1- 4H —)• 0, in the limit of a —)• r/-—r—, (4.7) 

z \ Ai Ap 

i.e., the effective mass of the fluid becomes equal to its wavenumber squared in the limit of the 
background geometry becoming de Sitter. Thus, when considering the issue of the decoupling 
of the fluid modes at the second transition, one must consider the full fluid effective mass H 22 
instead. As we will see in the next subsection, we shall be primarily interested in the scalar 
modes with wavenumbers k > femin- Using the value of fcmin obtained in Sec. 2.2, we have 

r^22 > 4in ~ (‘S’oAo)^'^^ ~ > Hi2 + -X '12 ~ X!i2 ~ 4'^^’ 

where we used k^iin ~ (TqAo)^/^ and Tq ~ e^/Ao consistent with our choice of Ap = 1 at 
the second transition (also note that we are assuming Aq ~ Ai for simplicity). Therefore, 
it appears that the fluid and inflaton perturbation modes of primary interest will eventually 
decouple again early on during the second transition. 
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4.1 Adiabaticity of the scalar perturbations ; w = ^ 

The analysis done in Sec. 4 indicates that the fluid and inflaton modes dynamically decouple 
at early times, reducing to two uncoupled harmonic oscillators with time dependent frequen¬ 
cies. We now show that if the fluid is radiation (w = 1/3), both the fluid and inflaton modes 
are in an adiabatic regime at early times, that is the effective frequencies -v/Oj/ of those modes 
evolve in timescales much longer than the frequencies themselves. This concept of adiabatic¬ 
ity replaces and generalizes the usual concept of modes being "inside the horizon". Indeed 
for both the fluid and inflaton modes, the effect of the background geometry on subhorizon 
perturbations diminishes as with ry —)• 0. This can be seen in Table. 1, where we showed 


O 22 = cik^ - — -4A = cie + 


This indicates 


[4z^A - {zB)'] = e + 0{rf). 

(4.9) 


^22 


u/i 

^^22 


“11 


< 1 


(4.10) 


as ry —7- 0. Therefore in this limit the scalar modes behave like simple, time-independent 
harmonic oscillator for which the general state furnishes a natural choice of initial state, a 
generalization of the usual Bunch-Davies vacuum. 

Although all modes become adiabatic sufficiently early in time, not all modes will have 
become adiabatic by the time the Planck scale is reached, which is the earliest time our model 
can be trusted. We look into this matter by computing the adiabaticity ratio (4.10) evaluated 
for /cmin when the fluid becomes Planckian. Using the background perturbative solutions 
found in Sec. 2.3, we And (for w = 1/3) 




11 




3/2 

11 


{kn 


a; 


3/4 


Q.', 


22 




3/2 

22 


{kr. 


10 A, 


1/4 


(4.11) 


where Aq = U((^o)- Note that here we have specialized to the class of solutions for the inflaton 
field for which the held freezes at early times. Since the inflaton field is expected to be at the 
GUT scale during inflation, we And 




11 


54 


3/2 

11 


(fc) 


> 10 


-6 




22 


54 


3/2 

22 


(fc) 


> 0.1, for all k<k^ 


(4.12) 


Evidently, the fluid modes with wavenumbers smaller than k^\^ will not enjoy a strong adi¬ 
abaticity at the Planck scale. We will therefore be interested in having the observable scales 
correspond to wavenumbers larger than kadn- 


4.2 Adiabaticity of the tensor perturbations ; w = ^ 

Just as in the scalar perturbation case, the tensor modes are in an adiabatic regime at early 
times. For tensor perturbations to be in an adiabatic regime at early times, the frequency of 
the gravitational waves ui defined in Eq. (3.33) must satisfy 

to' 

^ <1. (4.13) 

Using the background solution for the scale factor found in Sec. 2.3 one finds 

Lo‘^{k) = k^ — — = k^ + (4.14) 

a 
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indicating that uj‘^{k) —)• k'^ for all tensor modes as ?? —>■ 0. However, we mnst see for what 
modes the conditions given in Eq. (4.13) are met by the time the energy density of the flnid 
becomes Planckian. A short calcnlation similar to the one done for the scalar modes shows 


W (^min) 
^^(/umin) 


A, 


1/4 
0 > 


(4.15) 


where the expression is calcnlated when the flnid is at the Planck scale. Therefore for all 
gravitational waves with wavennmbers k > femin, there exists a natnral choice of gronnd state 
at the initial time when the flnid energy density is Planckian. 


4.3 Non-adiabaticity of perturbations ; tc / 4 

In the case of a flnid with an eqnation of state parameter w ^ 1/3, a short calcnlation nsing 
the resnlts of Sec. 2.3 gives 


= c^k"^ — 


2 — 6w 


(1 + 3rc)2 

2 ; 2 
OJ = k — 


7 ? ^ + Oikorf), 
2 — 6w 


Hii = k — 


2 — 6w 
(1 + 3u))2 


7] ^ + O{kori^), 


(1 + 3tt;)2 


T] + 0{Aor]^). 


(4.16) 


Here one cannot achieve initial adiabaticity in the sense of (4.10) for all modes of interest, 
particnlarly for modes near /Cmin which are of primary interest in this paper. 

The non-adiabaticity for a mode is an indication of its interaction with the backgronnd 
geometry. In snch cases one no longer has a natnral choice of gronnd state. One conld take 
different approaches to parameterize snch initial states as excited states, the core idea being 
to constrain snch states nsing backreaction considerations of the energy density stored in 
these initial qnanta. Computations nsing excited initial states are more specnlative dne to 
onr ignorance of the correct initial occnpation nnmbers for pertnrbation qnanta. Therefore 
we refrain from constrnoting nnmerical examples for snch scenarios in sections 6.1 and 7.1. 


5 Curvature and entropy perturbations in this model 


In the cosmological theory of small flnctnations, conservation laws can exist for certain gange 
invariant variables in some circnmstances. The existence of conserved qnantities is the primary 
reason why predictions of some inflationary models can be tested despite onr ignorance of 
what went on immediately after inflation. Of primary interest in this regard is a gange 
invariant fnnction known as the comoving cnrvatnre pertnrbation TZ dehned by 


n = ^-n 


Sq _ I 

a{p^ + pO) a ’ 


(5.1) 


where 5q and 5u are the three-momentnm and the three-velocity potentials associated with a 
flnid [see appendix A]. The comoving cnrvatnre pertnrbation gives the initial conditions for 
the temperatnre and polarization flnctnations that we observe today in the CMB [5]. Taking 
TZ to the Fonrier space nsing TZ{x, r/) = f d^k/{2Tr)^^‘^ TZ{k, one can show [7] 


T-IrP' f k"^ 

TZ\k,v) = -3^S{k,v) + 0(:^ 


P 


(5.2) 
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where S{k,rj) is the Fourier transform of a gauge invariant function called the entropy per¬ 
turbation, 


S = -H 


5p 


pO/j’ 


(5.3) 


where 5p and 5p are the pressure and energy density perturbations. This function is the gauge 
invariant measure of the non-adiabaticity of scalar perturbations. What we are primarily 
interested in is whether the comoving curvature perturbation for the fluid-inflaton model is 
conserved on super horizon scales. 

It is a remarkable fact that for any background cosmology, one can hnd solutions of 
the perturbation equations (3.27) for which TZ(k) is conserved as long as it remains in the 
superhorizon regime, i.e. as long as k Ti [4, 30-32]. This fact remains valid regardless of 
the specihcs of all subsequent cosmological eras. These perturbations, known as the adiabatic 
perturbations, can be intuitively realized as local spatial curvature fluctuations inducing a 
time-delay effect in the background geometry [33]. All linear fluctuations in single component 
cosmological models are known to be of this kind. On the other hand, for more complicated 
cosmological models, such as multiheld inflationary models, non-adiabatic solutions can exist. 
For these solutions, S{k) will generally not vanish in the superhorizon regime. Non-vanishing 
of S{k) causes the comoving curvature perturbation TZ{k) to change even in the superhorizon 
regime. Without a quantity that remains conserved in the superhorizon regime, predictions 
of an inflationary model are only testable if one chooses a specihc theory of reheating. 

Luckily, the entropy perturbations are unimportant in our model. We establish this fact 
by showing that the entropy perturbation function S reduces to its single held counterpart 
after the second transition. Using the dehnition (5.3) we have 


5 


n 


^Pf T ^Pip 


6pf + 5p^ - 

pO/ 


(5.4) 


where ^Pf{p) and 5pf[p) are the energy density and pressure perturbations of the huid (in- 
haton), while p^ and p^ are the overall background energy density and pressure. A short 
computation using the Friedmann equations (2.8) and (2.9) gives 


P = 



/ 

p = 



e(3-2r?). 


(5.5) 


where e = and p = —H/{2HH) are the two main slow-roll functions parameterizing 

the degree to which the background de Sitter symmetry is broken. It is also convenient to 
dehne 


2^3 

(5.6) 

where the subscript (p{f) on derivatives of the Hubble parameter denotes the contribution 
due to the inhaton (huid). 

Next we examine the huid and inhaton perturbation functions. The inhaton perturbation 
functions are [7] 

6p^ = dX + V^dip, 6p^ = 6X -V^5ip, (5.7) 


0 __^ 

~ 2iL3 ’ 




Sin — 
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where we defined 6X = ip'S^p' ja? — cf)ja?. The fluid perturbation functions can be derived 

from the linearized Einstein equations [8] 

= 5pf = cl5pf, (5.8) 

which using the constraint equation (3.14) reduces to 


5pf = - 


a^c‘1 L 


+ 0^13(1) 


(5.9) 


Using Eqs. (5.6) to (5.9), we write 5 as a power series in the fluid slow-roll parameters rjj 
and Ef, 

5 = + ... (5.10) 


with 


and 


c(o) _ 1 j,/ 

-^3(3-27/0)'^ 


(5.11) 


5^^) = — 

n 


+("^99 ~ ^/) 


2« - 3) 
3(3 - 27?0)e,, 

1 2 
^ ~ 3 




2r?: 


!(/_£/ 

(3 - 2r/0)2 \p0 




H 


3 - 




Tl! _ 

^ n 






1 + 


3ci 


3-277O 


n'f 


(5.12) 


Here we have used the dehnitions 




and 


nf = - 


V2z 


(5.13) 


which are the comoving curvature perturbations in absence of the other species. Judging from 
the slow-roll parameters appearing in Eq. (5.12), one can see that after the second transition 
when <C and Ef <C S converges to which is the entropy perturbation function 
of a single inflaton held. This suggests that as the huid energy density becomes subdominant 
to that of inhaton, the entropy between the two species diminishes and the overall entropy 
perturbation function reduces to its value in the absence of huid. This indicates that our 
model does not produce non-negligible entropy perturbations, just as in the case of single 
held inhation. 


6 Power spectrum of the scalar perturbations 


The power spectrum of the comoving curvature perturbation is the primary tool by which 
predictions of inhationary models are experimentally tested. Using the dehnition (5.1) we 
can write 


7 ^ = '0 — “77 


6qf + 6q^ 

+ 2/3 


ip 


n 

+ 2/3 


/ r I 

p dp + Cs - 0 

a 


_ 

(^'2 + 2/3 p >‘^ + 2 p ' 


^'21 


( 6 . 1 ) 


where the sum of momentum perturbations Sqp + dq^p was obtained from the right hand side 
of Eq. (3.16), being —Tfl2. Also the gauge invariant variables v and vr were introduced in 
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Eq. (3.18). Note that 7r/(^ and vl\/2z would be the comoving curvature perturbations in 
inflaton or fluid-only models. 

Since we are interested in understanding how primordial quantum fluctuations seed the 
formation of CMB anisotropies, we need to compute the quantum correlation functions of 
the comoving curvature perturbation TZ. The presence of derivative interactions between v 
and TT in action (3.19) renders the quantization procedure non-trivial. The details of this 
procedure is provided in appendix C. It turns out that one can write the quantized modes of 
the comoving curvature perturbation TZ as 


7^k(7y) = W^(r/)nk(r/). 


Here 


W = 


A-i 

w ' 2 + 2/38 




\v5'2+2/3' 

and the quantized modes for the generalized solutions IIk(??) are defined as 


( 6 . 2 ) 

(6.3) 


flkiv) = ^A{k, ri)aA,\, + n^(fc, ??)«A _k> (6-4) 

where similar to the notation used in Sec. 4 





(6.5) 


with the set {111,112,11)^,112} of four linearly independent solutions spanning the four di¬ 
mensional space of solutions for the coupled system of differential equations (3.27). Also, 
the annihilation operators a^k eliminate the quantum mechanical ground states |H)/i for the 
inflaton and the fluid perturbations. The existence of natural ground states is guaranteed by 
virtue of the adiabaticity of the scalar modes at early times discussed in Sec. 4 for the case 
of inflaton coupled to radiation, which is our primary interest. 

Given the fact that there is a natural choice of ground state for the inflaton and fluid 
perturbations, we compute the power spectrum of the comoving curvature perturbation in 
the initial state 

|X) = |Hi)® IH 2 ), (6.6) 

where |Hi(2)) is the initial ground state for inflaton (fluid) perturbations. The power spectrum 
of the comoving curvature perturbation then becomes 

{I\TZl{rjout)nk'{r]out)m = + k')7^7^(fc) (6.7) 


where ?7out is conventionally taken to be the conformal time for which k = Ti, and 


= Y, [w^UAik]) (w^U*Aik)) 

A 


E 


A L 






(y'2 + 2/3 


TTA{k)\‘^ f 2/3 y |t/A(fc)p ^ f 2(3ip'‘^ TTA{k)v*^{k) ] 

e ^W^ + 2pJ 2z2 (^'2+ 2/3)2 JJ’ 

( 6 . 8 ) 


where all the background and perturbation functions are evaluated at ryout- Note that the 
scalar power spectrum receives most of its contribution from the first term in the square 
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bracket, since —)• 0 with the fluid energy density becoming subdominant to that of 

inflaton. 

It is often useful to work with a dimensionless power spectrum defined as 


A 


9 T\j 




(6.9) 


Since the standard picture of inflation predicts a nearly scale invariant power spectrum, A^ 
is expressed using a simple ansatz 


A 


2 

'll 


= A, 


k 

h) 


( 6 . 10 ) 


where 


. 1 
As = —^ — 
" 87r2 e 


is the value of the power spectrum evaluated at a pivot scale /c* and 


( 6 . 11 ) 


Us — 1 = 2r/ — 4e + terms higher order in slow-roll, (6.12) 

is a measure of the scale dependence of the scalar power spectrum [7]. 


6.1 Case study 

We now numerically evaluate the power spectrum (6.9) for the case of radiation coupled to 
an inflaton with a potential. As discussed previously, specializing to radiation has the 
advantage of bringing about the adiabaticity of perturbation modes at initial times. 

As already discussed in Sec. 2, we shall fix the initial energy scale of the radiation at 
the Planck scale and specialize to the class of background cosmology for which ip approaches 
a constant value in the —?■ 0 limit. In this case, the remaining free parameters are the 
initial value of the inflaton field and the coupling parameter of the inflaton potential m?/2. A 
super-Planckian field variation of order Api ~ 15 during inflation together with the measured 
values of A^ and its scale dependence parameter Ug require 10“® ^ ^ 10“^ [7]. Knowing 

the range of these two parameters, we look for values that result in a suitable fit to the 
reconstructed values of the amplitude Ag and the scale dependence parameter Ug for the p>‘^ 
potential model at the pivot scale = 0.05 Mpc“^, as reported by the Planck collaboration 
[ 10 ]. 

In Fig. 6 we produce a plot of A^ for a range of low multipoles 2 < Z < 28, assuming 
that inflation is followed by instantaneous reheating. Different plots of the scalar power 
spectrum correspond to different values of an observationally relevant parameter 91, which is 
the dimensionless ratio of the fluid to inflaton energy densities when the dipole scale I = 2 
exists the horizon. Fixing the value of m, we change 91 by changing the initial values of the 
inflaton field. 

A number remarks are in order here: 


• In Sec. 4.1 we argued that modes with wavenumbers k < /Cmin have never been in an 
adiabatic regime since the Planck era. This means that the dipole scale mode is not 
initially in an adiabatic regime for 91 > 1. For the case of 91 = 10, one can show that 
multipoles I > 5 are initially in an adiabatic regime. 
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Figure 6: The power spectrum versus k for radiation coupled to an inflaton with the potential V(tp) = 
12 for 2 < I < 28 assuming instantaneous reheating. The parameter 91 is the ratio of the fluid energy 
density to the inflaton energy density when the 1 = 2 mode left the horizon. 


• The coupling of inflaton to radiation produces two main features in the scalar power 
spectrum: large amplitude suppression followed by small oscillations for low multipoles. 
The oscillations die off rapidly as k increases and the scalar power spectrum quickly 
converges to the standard single held result. The period of oscillations, which is larger 
for larger values of IH, can extend over a few multipoles. Nonetheless, the amplitude of 
oscillations does not exceed a few percent. 

• Observational data for I < 5 are not reliable due to a high value of cosmic variance. 
For values of IK > 1, one can expect deviations from the standard single held result for 
Z > 5, as indicated in Fig. 6 for IH = 10. Even though the small oscillations produced 
in this model do not appear to be detectable, the sharp suppression at the end tail of 
the power spectrum could very well be. Nonetheless, this suppression is only slightly 
outside of the range of multipoles with large cosmic variance. Also one should bear in 
mind that these modes are barely adiabatic at the initial time. 

• Oscillatory features for low multipoles have previously been discovered in a variety 
of models, including trans-Planckian models [34-37], axion monodromy models with 
additional instanton induced oscillatory features in the potential [38], models with brief 
non slow-roll phases [15, 39-42], models of cascade inhation [43, 44], and models with a 
sudden change in the speed of sound [45-47]. In our case, the oscillations are produced 
while the adiabaticity for all modes within the observable range as well as the slow-roll 
nature of the background dynamics are preserved. 

• Sharp suppressions of the scalar power spectrum for small multipoles have been pre¬ 
viously suggested in cutoff models [48], models of cascade inflation [43, 44], as well as 
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models with an initial non-slow-roll dynamics [15]. Suppressions in the former class of 
models extend over a larger range of multipoles than in our case. However, the ones 
observed in the latter class of models are quite similar to ours. 

• The features observed in this particular case study are in a good agreement with the ones 
seen in the previous analyses of this model. In particular, both the oscillatory features 
as well as the sharp suppressions pertaining to the low multipole scalar perturbation 
modes were noticed in [49] and [50]. The large suppressions were also noticed in [51]. 

• Finally we comment on the effects of a non-trivial reheating scenario on the afore¬ 
mentioned results. A non-trivial reheating scenario, such as the canonical reheating 
mentioned in Sec. 2.4, results in a smaller ratio Oj/ao of the initial scale factor and 
the scale factor today. Therefore, the dipole momentum scale today would correspond 
to a larger momentum scale when the fluid energy density is Planckian. Larger val¬ 
ues of k exit the horizon at later times, when the ratio of fluid energy density over 
the inflaton energy density is smaller. Hence lowering the temperature of reheating in 
non-trivial reheating models correspond to having smaller deviations from the standard 
power spectrum at low multipoles. 


7 Power spectrum of gravitational waves 


In this section we compute the power spectrum of the stochastic background of gravitational 
waves produced in this model. The quantization procedure for the gravitational waves in this 
model is quite standard and can be found in many references [5-8]. 

The power spectrum of the gravitational waves is defined as 


{J\kj{k, rjout)h'^{k', Vout)\J) = S^{k + k')Vh{k), 


(7.1) 


where r/out is the conformal time when k = 7i, and \J') is the Bunch Davies initial state. Here 
we have [5] 

|h(fc,r/out)P 


Vh{k) = 4i 


(7.2) 


a^illont) 

Just as in the case of scalar power spectrum, it is convenient to define a dimensionless 
tensor power spectrum A| using 

^ (7.3) 

In the standard inflationary picture, due to the near scale invariance of the spectrum of 
tensorial fluctuations, one normally expresses A| using a simple ansatz 


^i = At 


k\'^t 


where 


A, = 


TT^ 


(7.4) 


(7.5) 


is the value of the gravitational waves power spectrum evaluated at a pivot scale A;*, and 

nt = —2e -|- terms higher order in slow-roll, (7-6) 

is a measure of the scale dependence of the tensor power spectrum [7]. 
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7.1 Case study 

We numerically compute the tensor power spectrum from Eq. (7.3) for the same model 
studied in Sec. 6.1. We choose the Hubble parameter [12] 

ET = 3 X W — X 10"® « 4.24 X 10"^ (7.7) 

V 0.1 ^ ^ 

for I ~ 80, which corresponds to a scalar to tensor ratio of r = 0.2. Using this value of the 
Hubble parameter, we calculate the amplitude At of primordial gravitational waves to be 

At = ^ 3.64 X 10"^°, (7.8) 

near I ~ 80. As in Sec. 6.1, we select values of the inflaton mass m and the initial value of the 
inflaton field that give both an At close to the measured value near I = 80 and nt ~ — 2e, as 
required from the theoretical considerations. Fixing the value of m, we vary the initial value 
of the inflaton field to produce various values of 91. Fig. 7 presents our result for different 
values of 91. 
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Figure 7: The tensor power spectrum versus k for a model with radiation coupled to an inflaton with 
the potential V{(p) = 1/2 for 2 < Z < 28 assuming instantaneous reheating. The initial value of the 

inflaton field is slightly changed to produce different values of 91. 

We now make a number of remarks with regard to this result: 

• For tensor perturbations, one can show by explicit calculation that the dipole scale is 
initially in an adiabatic regime even for 91 = 10. Recall that for 91 > 1, the dipole scale 
is smaller than A:min- 

• The tensor power spectrum exhibits the same qualitative features as the scalar power 
spectrum; namely small oscillations and sharp suppression for order unity multipoles. 
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The amplitude of oscillations however seem somewhat larger in this case, though not 
large enough to merit a phenomenological signihcance. Both of these features become 
insignihcant for / > 5. 

• Similar oscillatory features were noticed in a study of trans-Planckian models in [52], 

• The observed features are in a good qualitative agreement with the ones found in an 
earlier study of this model in [49]. 

• As in the case of scalar spectrum, having a non-trivial model of reheating corresponds 
to less deviations from the standard result in the observable window of scales. 

8 Concluding Remarks 

We recently learned that the mixed radiation-inflaton cosmological model studied in this 
paper had been previously studied with regard to the issue of the robustness of inflationary 
predictions hrst in [49] and subsequently in [50] and [51]. In [49], the authors computed the 
primordial scalar and tensor power spectra, as well as the CMB angular TT and TE power 
spectra. They showed that both the primordial power spectra and the CMB angular power 
spectra are suppressed at low multipole moments. They also found both the scalar and tensor 
primordial power spectra to have oscillatory features at low multipoles. A similar analysis was 
subsequently performed in [50]. There the authors computed the scalar power spectrum as 
well as the CMB angular TT and TE power spectra. Their analysis conhrmed the conclusions 
of the previous analysis performed in [49]. Additionally, in [51] the primordial power spectrum 
of a test scalar held in a hxed radiation-deSitter background was computed, using which the 
CMB angular TT power spectrum was approximated. Both the primordial and the CMB 
power spectra were found to be suppressed at low multipoles, though no oscillatory features 
were noted. 

In this work, we expand the previous analyses by a more detailed study of the background 
dynamics, as well as a complete treatment of the scalar perturbations. We note here that in 
all of the previous studies the scalar perturbation due to the radiation huid was not properly 
included. As such, the possibility of the generation of the entropic perturbations was not 
addressed. In this work, by virtue of the scalar perturbation action derived in Sec. 3, we 
were able to identify the exact form of the linear order couplings between the radiation and 
inflaton perturbation modes. The study of these coupling functions together with the study 
of the background transition epochs were instrumental in proving several key properties of 
the perturbations. These include the arguments provided in Sec. 5 for the suppression of 
the entropy perturbation modes. Nonetheless, our analysis conhrms the validity of the major 
conclusions previously drawn in [49-51]. In particular, the results presented in Sec. 6.1 for 
a case study of the scalar perturbations is in a very good agreement with the conclusions of 
[49, 50], and to a lesser degree with those of [51]. Also, our results for a case study of the 
tensor perturbations presented in Sec. 7.1 are in good agreement with the ones found in [49]. 

Finally, we emphasize that our choice of a radiation fluid for the pre-inflationary phase is 
motivated by the fact that in this model the perturbations at sufficiently early times are adi¬ 
abatic and so a natural choice of initial vacuum state exists. However, just as with standard 
inflation, initial non-vacuum states may be physically relevant for a variety of reasons, includ¬ 
ing the effect of trans-Planckian physics or the effects of a phase that preceeds the adiabatic 
radiation-inffaton phase. For standard inffation, one typically parametrizes trans-Planckian 
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effects on inflationary dynamics by deviations from the initial Bunch-Davies vacuum state^ 
[see e.g. [53, 54] for some of the earliest approaches to this problem]. While this is an 
interesting issue, we have not addressed it here as it is beyond the scope of this paper. 
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A Review of stress-energy tensor perturbations 


In this section we review how a general stress-energy tensor for a fluid can be linearly per¬ 
turbed. For a perfect fluid, the stress-energy tensor can be written as 


Tab = PQab + {p+ p)UaUb, 


(A.l) 


where p, p, and are the pressure, energy density, and the four-velocity of the fluid. An 
everywhere isotropic Universe requires the components of this tensor to only depend on time. 
However, allowing small spatial dependence for the fluid functions, 

P{V, x) = + 6p{rj, x), p(r/, x) = p^{r]) + 6p{r], x), 

u^(?],x) = uS(rj) +du‘^(?],x), (A.2) 


one is left with a slightly inhomogeneous and anisotropic Universe described by the linearly 
perturbed Friedmann-Robertson-Walker metric [5] 


ds^ 


-a(7]f 


(1 -F 2(l))drf‘ 


2Bidx'^drj — [(1 — 2ip)5ij + 2Eij]dx^dx^ 


(A.3) 


where 4>, Bi, and Eij are perturbation functions all being small compared to unity. 

We then dehne the background fluid four-velocity vector Uq to be that of the frame in 
which the three-momentum density vanishes. Recalling that pabU^'u^ = — 1 to all orders in 
perturbation, we dehne 


uo = 



6u = 



(/> u* - R* 
a ’ a 


(A.4) 


where (u® — B^)la is the huid three-velocity vector to leading non-vanishing order for some 
three-vector v. 

With these dehnitions, the stress-energy tensor (A.l) becomes (to linear order in per¬ 
turbations) 

^00 = 0,‘^Ip^ + 6p + 2p®(/)], 

Fio = [p^{Bi - Vi) - p^Vi] , 

Tij = [p^dij +5p + 2p^{Eij - 'tjj)]. (A.5) 

Finally we dehne a useful additive function known as the three-momentum density perturba¬ 
tions hg® using 

hg® = (p° -|- p^)av'^. (A.6) 

®An initial thermal state could be a well motivated choice. See [19-21] where this choice of initial state is 
considered. 
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This quantity is proportional to the time-space component of the stress-energy tensor, 

(A.7) 

a 

For scalar perturbations, one can dehne a three-velocity potential as 5u^i = avi. This in turn 
gives rise to the three-momentum density potential 6q dehned as 

6q^i = (p° -F (A.8) 

B Review of gauge transformations and gauge invariant functions 

In this section we will review the basics of gauge transformation in cosmological perturbation 
theory. We will then provide the gauge transformation relations for the perturbation functions 
used in this paper. For an excellent discussion of this concept see [8] or a more recent review 
article [7]. 

Gauge transformation in cosmological perturbation theory is concerned with how per¬ 
turbation functions change under changes in coordinate systems dehned on some spacetime 
manifold. Given a set of coordinate system (r/, x*), we dehne an inhnitesimal shift in coordi¬ 
nates as 

fj = r] + X* = X*-h X*), (B.l) 

for some smooth functions and For scalar perturbations for some smooth 

function The line elements dehned in Eq. (3.6) are geometrical invariants, independent of 
a choice of coordinate system. We should therefore have 

ds'^ = gabiv, x')dx°-dx^ = gabifj, x*)^ix“dx^ (B.2) 

where gabiv, x^) and dx^" are the metric and differential 1-form dehned using a new coordinate 
system (fi,x^). We can use Eq. (B.2) to derive the changes in scalar metric perturbation 
functions, 

B = B + 

E = E-C 

(B.3) 

Tensorial metric perturbations hij are dehned to be traceless and transverse, rendering them 
gauge invariant. Imposing the traceless-transverse condition ensures that no piece of hij 
transforms as a scalar or a vector under the coordinate transformation (B.l). 

In addition to the metric perturbation functions, scalars dehned on spacetime also un¬ 
dergo transformations upon changing coordinate systems. For a scalar function a(r/,x*), we 
can dehne its perturbation as 

Sa{g, X*) = a{r], x*) - ao{g), (B.4) 

where ao(h) is the background value of a{r],x^). Keeping a{g,x^) invariant under the coor¬ 
dinate transformation, we get 

6a = 6a — (B-5) 
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The scalar functions important in this paper are the fluid three-velocity potential 6, the 
inflaton perturbation 6ip, the three-momentum density potential 5q for the fluid and inflaton, 
and the energy density and pressure of both species. These scalar functions transform under 
an arbitrary coordinate transformation as 

e = e- 

5(p = 5(p — 

Sq = Sq - aip^ + , 

Sp = Sp-p'^e, 

Sp = 5p-p'°f. (B.6) 


Using the gauge transformation relations derived above for the metric perturbation func¬ 
tions and scalars, one can form gauge-invariant functions that are linear in these perturbations. 
The ones used in this paper are the mukhanov variables and tt for the fluid and inflaton, 


u=—{e- 2 zi!), TT = aSip + C4>, 


the comoving curvature perturbation TZ, 


n = ii;-n 


6q 


a{p^ +p^) ’ 


and the entropy perturbation function S 


^ p'u p/u 


(B.7) 


(B.8) 


(B.9) 


C Quantization scheme for the scalar perturbations 


In this section we show how to properly quantize the scalar perturbations of our fluid-inflaton 
model. Our starting point is the scalar perturbation action 


<52^ 


dr] / d^x - 


n'^n' - Vin^ctv^n + n^rn + n'^An 


where the matrices are defined as 


n = 

r = 



/ ^i + l[Az^A-{zB)'] 







(C.l) 


(C.2) 


We have suppressed dependences on p for all background functions above. The Euler-Lagrange 
equations of motion for this action are 


n" - [c^v^ + r - ^A'] n + An' = o. 


(C.3) 
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Next we compute the Hamiltonian for this system. First note that the canonical mo¬ 
mentum P in this model is 

p _ ^ ^ 

where L is the Lagrangian appearing in the action (C.l). Employing the Legendre transfor¬ 
mation, the Hamiltonian becomes 


Hscalar = J dx^ [P^H' - C] 

+ S7in^cls7m-p^xn + n^{^x^\-T)u , 

which in Fourier space reduces to 



(C.5) 


Hscalar = j dk^ ^ [PfcP-fc + UlllU.k - Pfc . (C.6) 

Here we define 

^ = c^k'^ - T + ^X^X (C.7) 

( Our convention for Fourier transforming functions is X(a;, ry) = f dk^/(27r)^^^ X(k, 

The Hamilton’s equations of motion for each mode are 

Pju =-= —JjIIfc — —APfc, 

n; = = Pfc - ^xn,. (C.8) 

The Hamiltonian (C.6) is the starting point for quantization. After promoting Pfc and 
life to quantum operators, we impose the following equal-time quantization relations, 

[HjfcjHjTc'] =0) [P*fc)Pjfc'] =0) [HifcjPjfc'] = iSij6^{k + k'). (C.9) 


Quantum operators such as P and fl are often conveniently expanded in terms of the 
creation and annihilation operators and the mode functions that are solutions to the equations 
of motion (C.3) or (C.8). Note that the space of solutions of (C.3) is four-dimensional. 
Therefore, we can pick any two linearly independent solutions Iln and together with 
their complex conjugates to span the space of (Lagrangian) solutions. This leads to the 
expansion of the quantum operator 11 as 

nfc(??) = nn(fc, r])dA,k + n^(^, (c.io) 

where d and are the creation and annihilation operators, and summation over the index A 
is assumed. 

It remains to give a proper definition for the creation and annihilation operators d and 
aC To do this, we define the inner product on the space of solutions of (C.3) as 

(U,V) =i[U+V'-U'+V + U+AV], (C.ll) 

where -|- denotes the Hermitian conjugate of a matrix. Using the equations of motion (C.3) 
one can show this inner product to be conserved, i.e. dri{U, V) = 0, for all solutions, making 
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it an appropriate definition of inner prodnct on the space of Lagrangian solntions. We then 
require the solutions n^(A:,r/) appearing in (C.IO) to be normalized in the sense of 


{UAik, rj), IlBik, T])) = 5 ab- 


Conventionally we pick 


Ui{k,r]) 



n2(A:,r?) 



(C.12) 


(C.13) 


These two solutions are linearly independent, normalized in the sense of (C.12), and together 
with their complex conjugates form a basis for the space of solutions of (C.3) at early times, 
when the fluid and inflaton scalar perturbations are both decoupled and adiabatic. These 
solutions correspond to the fluid and inflaton perturbation modes being in some vacuum state 
at early times [8] . 

Finally we define the creation and annihilation operators as 


aA,k = {nA{k,r]),nk{r])), 

-{^A{k,v),tl^k{v))- 

(C.14) 

It is easy to check that 



[aA,k,a^B,k''^ = ^AsS^ik - k'), 

\^A,ki^B,k'\ 

(C.15) 

The operators aAk annihilate the ground states defined as 



= l^Ak), 

(C.16) 


k 


where lO^fc) denotes the vacuum state corresponding to the mode with momentum k. Ex¬ 
istence of a natural choice of vacuum at early times is the consequence of the adiabaticity 
of fluid and inflaton perturbation modes discussed in Sec. 4. Indeed the vacuum |fli(2)) is 
the analog of the Bunch-Davies state for the inflaton (fluid) perturbations, corresponding to 
mode functions with initial conditions (C.13). The vacuum state for the scalar perturbations 
is then defined as 

|X) = (g)|F!^). (C.17) 

A 
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